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We apply generic order parameter equations for the emergence of retinotopy between manifolds of different 
geometry to one- and two-dimensional Euclidean and spherical manifolds. To this end we elaborate both a linear 
and a nonlinear synergetic analysis which results in order parameter equations for the dynamics of connection 
weights between two cell sheets. Our results for strings are analogous to those for discrete linear chains obtained 
previously by Haussler and von der Malsburg. The case of planes turns out to be more involved as the two 
dimensions do not decouple in a trivial way. However, superimposing two modes under suitable conditions 
provides a state with a pronounced retinotopic character. In the case of spherical manifolds we show that 
the order parameter equations provide stable stationary solutions which correspond to retinotopic modes. A 
further analysis of higher modes furnishes proof that our model describes the emergence of a perfect one-to-one 
retinotopy between two spheres. 
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I. INTRODUCTION 

In a preceding paper JTl] we have analyzed a general model for the formation of retinotopic projections which is independent 
of geometry and dimension. In this paper we present applications of the general model, viz Euclidean and spherical geometries 
in one and two dimensions. To put these investigations into perspective we briefly recall their physiological motivations. But 
let us stress at the outset that our primary objective is not the biological modelling of retinotopic projections but the systematic 
analysis of a particular model thereof from a nonlinear dynamics point of view. In the course of ontogenesis of vertebrate 
animals well-ordered neural connections are established between retina and tectum, a part of the brain which plays an important 
role in processing optical information. At an initial stage of ontogenesis, the ganglion cells of the retina have random synaptic 
contacts with the tectum. In the adult animal, however, neighbouring retinal cells project onto neighbouring cells of the tectum 
[|2J. A detailed analytical treatment by Haussler and von der Malsburg was able to describe the generation of such retinotopic 
states from an undifferentiated initial state as a self-organization process In that work retina and tectum were treated as 
one-dimensional discrete cell arrays. The dynamics of the connection weights between retina and tectum was assumed to 
be governed by the Haussler-von der Malsburg equations which are based on modelling the interplay between cooperative 
and competitive interactions of the individual synaptic contacts. The nonlinear analysis was performed using the methods of 
synergetics, which provides effective analytical methods to study self-organization processes in complex systems 10 Ht]. 

Obviously, the description of cell sheets as linear chains with the same number of cells is an inadequate approach to the 
real biological situation. In a preceding paper we generalized the underlying Haussler-von der Malsburg equations to contin- 
uous manifolds of arbitrary geometry and dimension [1]. We performed an extensive synergetic analysis of these generalized 
Haussler-von der Malsburg equations. The resulting generic order parameter equations represented a central new result, and 
can now serve as a starting point to analyze in detail the self-organized emergence of one-to-one mappings in cell aiTays of 
different geometries. A short review of our generalization of the Haussler-von der Malsburg equations and the results of the 
corresponding synergetic analysis is provided in Sec.|II] In the subsequent two Sections we focus on one- and two-dimensional 
Euclidean manifolds. We show in Sec. |III]that the treatment of strings yields results which are analogous to those obtained for 
discrete linear chains in Ref. jst], i.e. our model includes the special case discussed by Haussler and von der Malsburg. However, 
our synergetic analysis is more general. Instead of discrete cell arrays with the same number of cells, we consider continuously 
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FIG. 1: Retina and tectum are represented as manifolds Mr and Mt, respectively, which are connected by positive connection weights 
w{t, r). The connectivity within each manifold is represented by cooperativity functions c_R(r, r') and crit, t'). 



distributed cells on strings of different lengths |@]. Furthermore, we do not restrict our investigations to monotonically decreas- 
ing cooperativity functions of strings. We investigate under what circumstances non-retinotopic modes become unstable and 
destroy retinotopic order. We show in Sec. |IV] that our generic order parameter equations also provide a suitable framework to 
describe the emergence of retinotopy between planes. For a certain superposition of two modes we demonstrate that taking into 
account the contribution of the higher modes leads to a sharpening of the retinotopic character of the projection between the 
two cell sheets. Finally, we analyse in Sec. |V]the formation of retinotopic projections for the biologically relevant situation of 
spherical geometries, i.e. manifolds with positive constant curvature. It turns out that the case of spheres exhibits remarkable 
similarities with the analysis of strings, especially regarding the generation of 1-1-retinotopic projections. 



II. THE GENERAL MODEL 



To make the present paper self-contained, we briefly review the essential results of our general model for the self-organized 
emergence of retinotopic projections between manifolds of different geometry in Ref. [1]. The two cell sheets, retina and tectum, 
are represented by general manifolds Ain and Mt, respectively. Every ordered pair {t, r) with t E A4t, r E Mris connected 
by a connection weight w{t, r) as is illustrated in Fig.[T] The equations of evolution of these connection weights are assumed to 
be given by a generalization of the Haussler-von der Malsburg equations 

w(t,r)^f{t,r,w)~^^ Jdt'f{t',r,w)~^^ jdr'f{ty,w), (1) 

where the first term on the right-hand side describes cooperative synaptic growth processes and the other terms stand for cor- 
responding competitive growth processes. Here Mt , denote the magnitudes of the manifolds, the total growth rates are 
defined by 

/(,..., = o + .(,,.,/di'/*V(M')c„(..V(,',0. (2) 

and a is the global growth rate of new synapses onto the tectum which represents the control parameter of our system. The 
cooperativity functions CT{t,t'), CT({r^r') represent the neural connectivity within each manifold. We assume that they are 
positive, symmetric with respect to their arguments and normalized. The cooperation strength depends on the distance between 
two points of the manifold. This requires a measure of distance, i.e. metrics, which in turn define Laplace-Beltrami operators on 
the manifolds. Their eigenvalue problems yield a complete orthonormal system ipx^ (t) , V'Ar {r), and the generalized Haussler- 
von der Malsburg equations are most conveniently transformed to this new basis. For example, the cooperativity functions are 
expanded in terms of these functions as follows: 

At Af£ 

The initial state of ontogenesis with randomly distributed synaptic contacts is described by the stationary uniform solution of 
the generalized Haussler-von der Malsburg equations wo{t,r) ~ 1. Its stability is analyzed by linearizing the Haussler-von 
der Malsburg equations ([1]) with respect to the deviation v{t, r) = w{t, r) — wo{t, r). The resulting linearized equations read 
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v{t, r) — L{t, r, v) with the Unear operator 

L{t,r,v) = -av{t,r)+ f dt' j dr' CT{t,t') CR{r,r')v{t' ,r') 
1 



2Mt 
1 

2Mr 



dt' 
dr' 



v{t',r)+ dt" dT"cTit',t")cR{r,r")v{t",r") 



v{t,r') 



dt'' 



dr"cTit,t")cR{r',r")v{t",r") 



(4) 



The eigenvalue problem of the linear operator (|4]l is solved by the eigenfunctions 

VAt A„ {t, r) = V'Ar (OV'Ah {t) (5) 

and the following spectrum of eigenvalues: 

~a — 1 \t = A/f = 

-« + \ ifL ft - 1) At = 0, ^ 0; Afl, = 0, At (6) 
~ct+fT„f^„ otherwise. 
The eigenvalue with the largest real part is given by Amax ~ —a + f^u /^t , where A^, A^ denote all those eigenvalues which 



A 



ArAfi 



FT fR. 

could become unstable simultaneously. Thus, the instability takes place when the global growth rate reaches its critical value 

«c = Re(/J,/i|). 



The linear stability analysis motivates to treat the nonlinear Haussler-von der Malsburg equations ([T]l near the instability by 
decomposing the deviation v{t,r) = w{t,r) — wo{t,r) in unstable and stable contributions according to v{t,r) — U{t,r) + 
S{t, r). With Einstein's sum convention we have for the unstable modes 

C/(i,r)-C/A.A5iV'A^(i)V'A5^(r), (7) 

and, correspondingly, 

S{t,r) = Sx^x^ipx^{t)iljx^{r) (8) 

represents the contribution of the stable modes. Note that the summation in (|8]i is performed over all parameters ( At ; A_r) except 
for (A^; A^), i.e. from now on the parameters (At; Xr) stand for the stable modes alone. With the help of the slaving principle 
of synergetics the original high-dimensional system can be reduced to a low-dimensional one which only contains the unstable 
amplitudes. The general form of the resulting order parameter equations is independent of the geometry of the problem and 
reads 



Uxi^xi - Aa^a^^ C/a^A5^ + ^a|;a|'a|" C^a^'a^^' C/aj"a- 

+ ^A|;A|'A|"A|"'t^A^'Ai^' C/aj"A-" Ux^'"X^^- 



(9) 



It contains, as is typical, a linear, a quadratic, and a cubic term of the order parameters. The corresponding coefficients can be 
expressed in terms of the expansion coefficients /a^, fxa of the cooperativity functions (O and integrals over products of the 
eigenfunctions ijjxT (t)^ i^Xa [f )'- 

(10) 
(11) 



-^A(i)A(2)...A(") = dx'>j/y^{x)il)x(i){x)il)x'.^){x) ■■■ ■4^xi-^){x), 



A(i)A(2)...A(") 



dx IpxM (x) V'aC^) (x) ■■■ '0A(") (x) ■ 



The quadratic coefficients read 



.Ay,A^ AJ „ „ tA^ 



(12) 



whereas the cubic coefficients are 



^r,^r'^r"^r"' 



1 



fx ^ ' " y A V 



j^j^ I\^> ^A^A^' Ja^"A^ 



R R R R 



I\u /\u //\u iii5xuXui 



1 



Jaj."aj."' ] + <^ [./at f\R + /a^' hi' ] IxI'Xt -^Af 'A„ ^ 2 



<5ArO '5a5,AJ.' + f\R)Ixi'\R 



+ '^A«0 5x^X1' (1 + Ixt) hl'X-j 



JJ'^T ^R 

A^''A^",A^'"A^'" 



(13) 
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As is common in synergetics, the cubic coefficients ( fTSl ) consist in general of two parts, one stemming from the order parameters 
themselves and the other representing the influence of the center manifold H on the order parameter dynamics according to 



^AtAh - ^AfAf^A^'A^i' t^A^AS^t/A^'A^i 

Here the center manifold coefficients H^I^^S xu,xu, are defined by 



(14) 



Aa" A" + Aa"'A" ' — A 



j\t jAh 

A^ A^ ' A^A^' 



1 



/Ml 



1 



(15) 



The order parameter equations (|9]l for the generalized Haussler-von der Malsburg equations ([T]) can now serve as a starting point 
for analysing the self-organized formation of retinotopic projections between manifolds of different geometry. 



III. STRINGS 



In this section we specialize the generic order parameter equations (|9]l to one-dimensional Euclidean manifolds of strings with 
different lengths Lt and Lr. We start with introducing the eigenfunctions in Subsec. UlI Al In Subsec. UlIBI we observe that 
the quadratic term vanishes and derive selection rules for the appearance of cubic terms. In this way we essentially simplify the 
calculation of order parameter equations as compared with Ref. Jsj]. Furthermore, we show that the order parameter equations 
represent a potential dynamics, and determine the underlying potential in Subsec. IIIICI A subsequent transformation from 
complex to real order parameters in Subsec. IIIIDI leads to constant phase-shift angles. Thus, in Subsec. UlIEI we reduce the order 
parameter dynamics to two variables which correspond to the amphtudes of two diagonal modes. These two modes compete with 
each other, until one of them vanishes. Within the potential picture this means that the stable uniform state becomes unstable 
and the system settles in one of the two potential minima, as is discussed in Subsec. IIII Gl After one of the diagonal modes has 
won, only such modes are excited which contribute to the sharpening of the diagonal. Approximately solving the Haussler-von 
der Malsburg equations in Subsec. UlI Hl leads to the following scenario: Above a critical global growth rate ac the uniform state 
W[){t,r) = 1 is stable. By decreasing the control parameter a, the projection gets sharper and sharper. Finally, if there is no 
global growth rate of new synapses any more, i.e. a ^ Q, the connection weights are given by Dirac's delta function. Thus, 
a perfect one-to-one retinotopic state is realized. We conclude this discussion of strings with comparing our results with the 
corresponding analysis of discrete linear chains in Subsec. HITll 



A. Eigenfunctions 

The magnitudes of the manifolds Mt and Mr are given by Mt — Lt and Mr — Lr, respectively. To avoid problems at 
the boundaries, we assume periodic boundary conditions, i.e. we consider retina and tectum to be rings with circumferences Lt 
and Lr, respectively. The eigenvalue problem of the Laplace-Beltrami operator for both manifolds reads 

ip\{x) ^ Xx'ipxix) , (16) 
with X — t,r, respectively. Using the boundary condition i)x{x) ~ ^x{x + L), this is solved by the eigenfunctions 

ipx{x) = ^exp (i^Ax ) , (17) 



where the eigenvalues are given by x\ ~ — 47r^A^/L^, with x E [0, L) and A = 0, ±1, ±2, .... Every eigenvalue x\, apart 
from the special case xo = 0, is two-fold degenerate. The eigenfunctions form a complete orthonormal system: 



L 



dxMx)rx'{x)=Sxx', MxWxix')=S{x-x'). (18) 



Q A— — oo 



Note that the orthonormality relation in ( fTSl l follows directly by inserting ( fTTI i, whereas the completeness relation is proven by 
taking into account the Poisson formula |7]. The cooperativity functions only depend on the distance, which is given by the 
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Euclidean distance \x — x'\, i.e. c{x, x') = c{x — x'). Their expansion in terms of the eigenfunctions ( fTTI ) corresponds to the 
Fourier series 



c{x -^') = J X! /^<5^P 



A=-c 



z — A [X — X ) 

1j 



(19) 



The expansion coefficients fx are independent of the sign of the parameters A, i.e. fx ~ as the cooperativity functions are 
symmetric with respect to their arguments: c{x — x') = c{x' — x). 

B. Synergetic Analysis 

To specialize the order parameter equations (|9|l to the case of strings, we have to determine the integrals (fTOl) and (fTTT i of 
products of eigenfunctions. With (fTTb we obtain 



-'a<i)A(2)... A(") 



A(i)A(2)... A(") 



1 

7l 
1 

7i 



<5a(i 



A(i)+A(2) + ...+A("),A J 



ri-2 



(^Ad 



A(i)+A(2) + ...+A("),0 ■ 



(20) 
(21) 



From these results one can immediately read off the special cases 



ly A" = <5a'+A",A , ^A' A" A'" — Y, ^>^'+^"+>^"'■^ ' J^' >" = ^A',-A" 



(22) 



In ( fT2b the first integral of (l22b occurs only with unstable values of A. As they can only differ by their sign, it follows A' + A" ^ 
A. Thus, we have x" = 0, i.e. the quadratic term (fT2] l in ^ vanishes: A^Tl'xT^,^!,, = 0. To arrive at a more concise 
representation, we split the cubic contribution in (|9]l in two terms according to 



^A|,A|'A|"A|"'f^A^'A]^' C/A^"A]i" Ux^"'X^^"' - Ki^xi^x^ + A'2,A^A5^ ■ 



(23) 



The first term Ki takes into account the contribution of the order parameters themselves, while the second term K2 represents 
the influence of the center manifold on the order parameter dynamics. Applying the integrals (|22] | leads to selection rules for 
the appearance of cubic terms. It turns out that only those sums lead to non-vanishing contributions where the sum of three 
unstable modes A" + A" + A" coincides with another unstable mode A" . Thus, both for retina and tectum only the following 
combinations are allowed: (A"', A"", A"'") = (A", A", -A") , (A", -A", A") , (-A", A", A") . With this selection rule the first 
cubic term is given by 



fT fR 

2LtLii 



^A^A^ ^A^ A^ ^ A^ A^ '^'^j^ "'"'^H ^""^^ ''^i^ 



A^ ■^i^ A^ '^j^ A^ "^j^ A^ -|-A^ -|-A^ lA^ 



(24) 



and the second cubic term reads 



X(5 



\rp \ A^ A^ /i 



x^'^x^" x^i^" {(/a^/4 + /a5./a|)'^a^'+a^,a 



(1 + fx0^xi+x'^^xi^yT^K^:x-^ + (1 + fl0^x-^+x'^,x-JyR^Kix-^ 



(25) 



The latter depends on the center manifold, which follows from ( fTSl l. and ( l22b to be 



H 



A- XI 



fT fB. 



AJi. A r 



-A^ '^A^ "^"'^H 



^A^ +A^ ,At ''^^ 



(26) 
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C. Complex Order Parameters 

We can therefore conclude that the order parameter equations for strings have the form 

J/AJAJi =/lAJA'A ([/,[/*) (27) 

with the complex function 

hx^x^iU, U*) = Axj^x^Uxj^xi + Axj^xiUl^^u^U-xi^-xi + Bx^x^Uxi^x^^U-x^xiUx^^xi ■ (28) 
Here we have introduced the coefficients 



iJA^A" - 



LtLr \ 2K\^^X'^ — A2A^,2A_^ 

7+ (72^^.0 - l)/2 7+ (70'2A?, - i)/2' 



LtLr 



4- 



2AAH.A^t — A2AH..n 2Kxi^x^^ — Ao,2X^ 



(29) 
(30) 



with the abbreviations 7'''^'^" := Jxt^Xr ^nd7 :— ^/^t'^r — fj^^f^^. Now we turn to the question whether the order parameter 
equations dZTl ) represent a potential dynamics. To this end we derived in Ref. f§\ a condition for the order parameter equations 
which allows one to conclude whether or not such a potential exists. The potential criterion reads 

dhx^x^^iU,U*) dK^',.'iU,U*) 

dU.^'.^. " dm... ' ^ ' 

T Ft T R 

which is, indeed, fulfilled for ( |28] l. Furthermore, we derived in Ref. [8] the following conditions for determining the underlying 
potential; 

. _ l dV{U,U*) 

^^^^ 2 5C/* ■ ^ ' 

Integrating (|32l ) yields the potential 



'^Bx}^x^^Ux^x'^JJ-x}^-x'fU-x}^x^^Ux}^-x^^ ■ (33) 



D. Real Order Parameters 

For technical purposes it has turned out to be useful to work with complex order parameters so far However, in order to 
investigate their contribution to a one-to-one mapping between the strings, we have to transform them to real variables. We 
construct at first the real modes from the eigenfunctions (|5]l, ( fTTl ) of the linear operator L according to 

cxT\R{t,r) = ^[vxTXR{t,r) +v^xT-XR{t,i')] ^ / cos ( -^At^ + -^A^r ) , (34) 

sxTXRit,r) = --[vxTXR{t,r) - v-XT-\R{t,r)] ^ ^== sinl —Xri + —XRr] . (35) 

Z ^/LtLr \-Lt J 

These two modes span a real subspace. If we set 

a = pcos'ip, b = psimp; p>0, "fp ^ {~tt,7t] , (36) 

the following relation results: 

/ 2tt 2tt \ 

acxj-XR{t, r) +bsxTX„{t,r) ^ pcosi —XTt+ —Xur ~ ip] . (37) 

\Lt l^R I 



(A^,A^,) = (1,1) 



(l.-l) 





FIG. 2: Diagonal modes of different orientation according to l l43t for the case (Ay, A^) = (1, 1) and (1, —1), respectively. Here the phase 
shifts are set to t/j = i/p = 0. 



Thus, the subspace consists of all phase-shifted functions of pcx^x^^ {t, r). Then the modes belonging to the unstable eigenvalue 
(A^, A^) are given by the modes c\u \^ {t, r) and c\u {t, r) as well as all phase-shifted functions. Rewriting the unstable part 



to real modes 



with real variables uj: 



U{t,r) ^ Ux}^\^^vx}^\i(t,r) + U-x^-xiV-\}^-\^^{t,r) 

, leads to 

U {t, r) = uicxj^x^ {t, r) + U2Saj, {t, r) + us,cxj^-xi {t, r) + U4Sx:^-xi (t, r) , 



f/Aj.A]^ = (ui-m2)/2, U-xj^-X'^^ {ui+iu2)/2, 
Ux^^xi^ = {u3 - iu4)/2 , U-xj^x^ = (w3 + m4)/2 . 
Inserting the transformations (|40] | into the complex potential (l33T l. we obtain the following real potential 



(38) 
(39) 

(40) 
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(41) 



The corresponding equations of evolution for the real order parameters are determined from (1411 1 according to 
—dV{ui)/duj . They read explicitly 



Ul 



U2 



U3 



U4 



Aa^; 



Aa^a]^ 



A A" A" + 



^A^A^^ 






Bxi^xi 


4 




4 


Ax^x-^ 








4 




4 


^A^A^lj 




huD- 


1 B^T^R 


4 


4 








Bx^x-^ 


4 




4 



{ul + ui) 

{uj + ul) 
(uf + ul) 
{ul+uj) 



Ul , 



U2 , 



U3 , 



U4 , 



(42) 



E. Constant Phase Shift Angles 

According to Eq. ( |37] | the unstable part ( [39] l can be written as a superposition of two diagonal modes of different orientation 



U{t, r) — ^cos 



77 COS 



27r 



27r 



— \rpt — A 



Lb, 



(43) 
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as is illustrated in Fig. |2] With ( |36] | and (IJTT i we have ui = ^cosV',W2 = C sin ■0 1^3 = r?cos</3,U4 = 77 sin Then the 
amplitudes of the phase-shift diagonal modes read 

i=^ul+ul, j^^^ul+ul (44) 

and the phase angles are given by tan ijj ^ ui/u2 , tan tp = u^^/u/^. From the order parameter equations (l42l i it follows ui/u2 — 
ui/u2 , U3/U4 — U3/U4 . Thus, performing a separation of variables and a subsequent integration leads to the relation ui/u2 — 
const , U3/U4 = const . Consequently, the four real equations ( |42] | are reduced to two equations for the mode amplitudes ^ and 

•q: 

V - (aa^aj^ + + V ■ (45) 

The corresponding potential is 

Vii, V) - -~^^ie + V') ^^{^ + V') - ^^ev' ■ (46) 

Thus, we have reduced the four complex order parameter equations dZTl i. (|28] | to two real order parameter equations i45[ with 
the potential (|46] |. 



F. Monotonous Cooperativity Functions 

So far our considerations are valid for arbitrary unstable modes (A^, A^) . According to the eigenvalue spectrum ^ the 
unstable modes are determined by the expansion coefficients f\ of the cooperativity functions. We therefore derive in this 
subsection some basic properties of these coefficients. In particular, we investigate the consequences of monotonically decreasing 
cooperativity functions for their expansion coefficients f\. As c{x) is positive and normalized, we conclude |/a| < 1. Using the 
Euler formula, the symmetry c{x) — c(—x), and integrating by parts, the expansion coefficients can be written in the form 

L/2 

/a = — ^ / c'(x)sinf^Ax) dx, (47) 



tta y ^ ' \L 



which makes the symmetry f\ = /_a manifest. If we assume monotonically decreasing cooperativity functions, i.e. dc/dx < 
for X G [0, L/2], we obtain /i > 0. Furthermore, we can show that /i is the largest expansion coefficient by considering the 
expression 







L/2 

'27r A 1 . /27r, 
—X 1 — - sm — Aa; 
Lj 1 A \ L 



dx . (48) 



Because of c' (a;) < it follows indeed /i — /a >0 VA^^Ojil. Together with |/a| < 1 the maximum eigenvalue of Q results 
to be A,„ax = —ci + f\Ji'- Hence in this case there are four unstable modes (A^, A^) — (±1, ±1), which corresponds to the 
result obtained in Ref. H. However, the most fundamental insight of our more general analysis is that the real order parameter 
equations (05]) are also valid in the case where the cooperativity functions are not monotonic so that any mode (Aj., A^) can 
become unstable. It is plausible that there is a pathological development in animals which corresponds to this case. 



G. Potential Properties 



We now analyze the properties of the potential ( l46b . For this purpose we restrict ourselves from now on to the unstable modes 
(A^, A^) = (±1, ±1), whose indices will be discarded for the sake of simplicity. Then the potential ( |46l ) reads 

v{^, v) = -^{e+ v') v') - ^ev' , (49) 
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FIG. 3: The potential V{(,, rj) according to Eq. l |49t with A > 0. The originally stable state ^ = 77 = becomes unstable. The system settles 
into one of the two minima, i.e. one of the two modes vanishes. The right plot shows the equipotential lines. Dark grey values correspond to 
small values of the potential V. 



where the coefficients A, B follow from d29] l. ( [30] l to be 

A = 

B = — 



7 


(2- 




LtLr 


\ 2A-A2,2y' 




7 


r 7+(7'^°-l)/2 7 


f (7°-^ - l)/2 


LtLr 


2A - A2,o 


2A - Ao,2 



(50) 
(51) 



From the condition W = we determine the extrema of V{^, rj) and assign them to a minimum, a maximum, or a saddle point. 
In the unstable region with A > the potential ?/), which is depicted in Fig. [3] has 

• a relative maximum at Pi (0, 0) , 



• two relative minima at ^2(0, ^-4 A/A) and P^{^J-AK/A, 0) . 

• a saddle point at Pi{^ --iK/ {A + B), ^-4A/(A + S)) . 



In the stable region with A < only the relative minimum ^ = ?/ = does exist. Initially, the system is in the stable uniform 
state WQ{t, r) = 1. This state becomes unstable if the control parameter a is decreased to the critical value ac — fi ■ The 
eigenvalue A,nax = ~a + fT fi becomes positive, and the minimum passes into a maximum. The system settles into one 
of the two equivalent minima, i.e. a symmetry breaking takes place. Thereby the two modes compete with each other and, 
subsequently, one of the two modes vanishes. Which of them vanishes depends on the initial conditions of ^ and If the 
condition 7^(0) > ^(0) is fulfilled, the ^-mode vanishes, and vice versa. 



H. One-To-One Retinotopy 



In the following we assume that, according to the potential dynamics discussed above, only one of the two modes remains. 
These two modes show a pronounced maximum for t ~ —r and t ^ r, respectively, as is shown in Fig.|2l To assess the influence 
of higher modes, we calculate the center manifold S{U) for the case ^ = and rj ^ and set U4 = without loss of generality. 
Then it follows from (l44l i that rj — U3, and we obtain for the unstable part ( [39l ) 



rr, . /27r 27r 

U [t, r) = rj cos 1 i 

\Lt Lr 

With the center manifold ( |26l l the stable part ([8]l, ( fT4b reads explicitly 



S{U) 



27 



VL^(2K^ 



47r ^ 

^2,2) \Lt 



A-K 



(52) 



(53) 
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(A^,A"„) = (2.2) 



(A^,A^) = (3.3) 





FIG. 4: Higher diagonal modes, whicli are excited by the unstable mode (A^, A^ 
of the (1, l)-mode. 



= (1, 1) of Fig.|2] They amplify the retinotopic character 



Thus, those modes are excited which strengthen the retinotopic character of the projection. With the help of the complex modes 
it can be seen that this is also the case for higher modes, i.e. for (Ati, A^) = (1, 1) exclusively the modes (2, 2), (3, 3) etc. are 
excited, which are depicted in Fig. |4] Therefore, we follow Ref. yl] and use an ansatz which contains only diagonal modes 
and insert it into the Haussler-von der Malsburg equations If we restrict ourselves to special cooperativity functions, the 
resulting recursion relations can be solved analytically by using the method of generating function. Note that our derivation of 
the solution of the recursion relations coiTesponds to the gravitating chain in Ref. 



1. Recursion Relations 



Motivated by the above remarks we investigate the Haussler-von der Malsburg equations for strings with the ansatz 

oo 

r) = VLtLr wxvx,-xit,r) , (54) 



w( 



A=-c 



where v\,^\{t, r) is defined by (O and (fTTI i. Thus, taking into account the decomposition (fT9] l of the cooperativity functions, 
the Haussler-von der Malsburg equations (HJ can be written as 



w{t,r) = ~a[w{t,r) - 1] + w{t,r)y^ LtLr ^ wxfl fxVx-xit^r) ~ p{w)w{t,r) , (55) 

A— — OO 

where we have introduced the abbreviation 

oo 

p{w)= J2 ^-i^ofjf^- (56) 



A=-c 



Inserting the ansatz (|54| | into ( fSSl l and comparing the coefficients of the Unearely independent functions vx,-x{t, r) yields 

Wo = - [a+p(w)] (wo - 1), (57) 

OO 

wx = -[a+p{w)]wx+ ^^-J^ifjff^ ^^0- (58) 

j=~oo 

As w{t, r) is positive lUl], we obtain thatp(u') > and a + p{w) > 0. Therefore, the stationary state is determined from ( |57] | to 
be Wo — 1. 



2. Special Cooperativity Functions 

We restrict our further considerations to the following form of the cooperativity functions (fT9] i: /o = 1 , /i , fj = 
for j 0, ±1. With the abbreviation 7 :— ff fi the previous result ( |56] | can be written as p{w) = 1 + 2jwiw^i, so that the 



11 



Q3 < Q2 < (Xi 




FIG. 5: Decreasing the control parameter a to smaller values, we read off from ( I64t and ( 165b that the connection weight converges to Dirac's 
delta function J66t. 



equations (l58ll for the stationary case reduce to the recursion relation 



(59) 



3. Generating Function 
To solve the recursion relation ( l59l l, we define the generating function 



(60) 



A=-c 



Multiplying ( |59] ) with + z ^ and performing the sum from A = 1 up to infinity yields to a linear algebraic equation which is 
solved by 

a 



E{z 



a + 2'ywl — 7W1 (z + z 
To determine the coefficients wx we expand the generating function (1611 1 into a Taylor series: 



£;(z) = - 



(a + 27k;^) w (zi — 



.A=l 



|Zl| < |z| < |Zl| 



with the abbreviations 



7W1 



Zl = ;^(l + V/T^4^) 



Comparing (|62] | with ( l60l ) by taking into account wo = 1 determines wi to be 



7 — a 



(61) 



(62) 



(63) 



(64) 



Thus, together with (|63T l it follows zi = wi, and the remaining coefficients turn out to be w\ = w['^', which is valid not only for 
A 7^ but also for A = due to wq — 1. 



By inserting the latter result into ( l54l i we obtain 

Wo(t,r) = 



4. Limiting Cases 



1 - wf 



1 - 2wi COS {^j^t ~ 22Lr] +w'i 



(65) 
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u'o(t, r) 



Bifurcation 



wo{t, r) = N5t^_r 
t 



wo{t, r) = NSt,r 

FIG. 6: Bifurcation in the vicinity of the instability point for the linear chain as analyzed in Ref. [J]. The quadratic arrangement of the two 
linear chains allows a concise representation of the connection weights. Dark gray means high connection weights between the corresponding 
cells t and r. At the uniform initial state all connection weights are equal. The bifurcation drives the system into one of the two possible states, 
which differ in their orientation. Decreasing the control parameter a to zero leads to a one-to-one retinotopy. Instead of a delta function in the 
continuous case, here the retinotopic order is described by Kronecker deltas. 



For wi = the stationary uniform state reduces to wqII, r) = IVt E [0, Lt) , r E [0, Lji). In the case wi — 1, i.e. wa = 1 VA, 
we find with the help of the Poisson formula tTJ 

wo{t,r)=5(^-^~^ . (66) 

Hence we have a situation which is illustrated in Fig.|5] If the control parameter a is in the neighborhood of 7, the connection 
weight is essentially uniform with a small maximum for t/ Lt — ^r/Lp. Further decreasing of a leads to a sharpening of the 
projection. In the case a ^ the projection becomes Dirac's delta function, i.e. a perfect one-to-one retinotopy is achieved. 
This means that the undifferentiated growth of new synaptic contacts comes to an end when the ordered projection between 
retina and tectum is fully developed. 



I. Comparison with Linear Chains 

Finally, we compare our results for strings with those of Ref. Js]] where retina and tectum were treated as linear chains 
consisting of N cells, respectively. In that reference, the order parameter equations read 



with the abbreviations 



U,j = [A - 7(2 - a)U,jU-,-j + (4 - 6' - 6")C/»-,[/-y](7y (67) 

7 + 7^-^ w 7 + (7^'°-l)/2 ,„ 7 + (7"'^-l)/2 
a = , h = , b = . (68) 

A22 ^1-20 A02 

The comparison of the coefficients A, B according to ( fSOb , ( BTT i with 7(2 — a), 7(4 — &' — b") exhibits two differences: the 
factor 1 / LtLr in A and B as well as the term 2A in the denominator. The absence of the corresponding factor 1 /N"^ in Ref. fsl 
stems from the circumstance that the eigenfunctions were not normalized there. Physically more interesting is the appearance 
of the term 2A in the denominator of A, B. The reason for this is that we have used the mathematically correct equation 
for determining the center manifold ( fTSl l according to Refs. whereas in Ref. \3\ the center manifold is adiabatically 

approximated by 5 = 0. However, this ad-hoc method for implementing the adiabatic approximation, which is frequently used 
in the literature, is only justified for real eigenvalues. As the eigenvalues of the strings are real, we deduce for the vicinity of 
the instability point the relation A = —a + 7 w 0. Thus, the coefficients ( fSOl ). ( BTI ) turn into those of Ref. |3] and the adiabatic 



13 



approximation 5* = can be applied here. 

Furthermore, our results for the continuous case are analogous to the results for discrete cell arrays. Also the transition to a 
perfect one-to-one retinotopy takes place in a corresponding way, as is illustrated in Fig.|6] Thus, we conclude that our geometry- 
independent model for the emergence of retinotopic projections developed in Ref. 1 1 ] contains as a special case the results of 
Ref. Jit]. In addition, we have extended the range of validity, i.e. the domain around the instability with A = 0, where the 
order parameter equations represent a quantitatively good approximation, as we have derived a more precise form of the center 
manifold (fTSl l. 



IV. PLANES 



In this section we extend our discussion to two dimensions where the cell sheets are assumed to be planes of side lengths 
, L2 and if-, L2, respectively. To obtain a consistent solution we assume again periodic boundary conditions, i.e. the cell 
sheets are modelled as surfaces of tori. We start with presenting the linear analysis in Subsecs. llV Al and II V B I Afterwards, in 
Subsecs. llVCI and lIVDI it turns out that we have to calculate in total sixteen order parameter equations where the quadratic 
term vanishes, as in the case of strings, and where again selection rules reduce the number of cubic terms. This order 
parameter dynamics turns out to be complicated as the two dimensions do not decouple in a trivial way. Therefore, we 
have to restrict our analytical discussion of the order parameter dynamics to physiologically interesting special cases. If we 
set all modes to zero except for one, we find retinotopy only in one dimension, which is shown in Subsec. IIV El In a next 
step we consider the superposition of two retinotopic modes and investigate the necessary conditions for their coexistence. 
Such a situation occurs, for instance, when the cooperativity function of the tectum is monotonically decreasing, whereas the 
cooperativity function of the retina is not monotonic. In Subsec. IIV Fl we show that taking into account the center manifold 
contribution or higher modes leads to a sharpening of the retinotopic character of the projection between planar retina and tectum. 



A. Eigenfunctions 



In the following we consider both retina and tectum to be planes with side lenghts Li and L2. The points on the plane 
are represented by Cartesian coordinates x — (xi, X2) , xi G [0, Li) , X2 G [0, L2)- The magnitude of the plane is given by 
M = L1-L2. The corresponding eigenvalue equation A'i/j{x) ~ xi'ix) is solved for periodic boundary conditions, i.e. ipj ixj ) = 
'i/'j [xj + Lj ) by the complete orthonormal system of eigenfunctions 



'^\{x) 



1 



: exp 



Xixi X2X2 



Li 



L2 



where A — (Ai, A2). The cooperativity function c{x — x') is expanded according to (O in this basis: 



c{x — x') 



LIL2 



A2)exp<^ 2m 



Al .Ao 



Ai(xi 



X2{X2 



Li 



(69) 



(70) 



Note that again the expansion coefficients /(Ai,A2) of the cooperativity functions are independent of the signs of the parameters 
Al, A2, as the cooperativity functions should be symmetric with respect to their arguments: c{x — x') = c{x' — x). This 
requirement and the linear independence of the exponential functions leads to /(Ai,A2) = /(±Ai,±A2)- From now on we assume 
that the cooperativity functions decouple with respect to the two dimensions: c{x — x') — ci{xi — x[)c2{x2 — ajj). As the 
individual cooperativity functions can be expanded according to 



J = 1,2, 



(71) 



the decoupling amounts to a factorization of the expansion coefficients: /( Ai , A2 ) — /ai /a2 ■ With this we allow for both isotropic 
and certain anisotropic cooperativity functions. This is an interesting feature as it is reasonable to assume that real cell sheets 
have a preferential direction. 



B. Instability Point 



We analyze which modes become unstable. According to (|6]l this depends on the expansion coefficients and the maximum 

■T fR 
A" ^A" 



eigenvalues are given by Aa" a" = —a + /Ju . By doing so we require that all unstable modes become unstable simultane- 



14 



ously. This requirement is due to the fact that the order parameter equations should be approximately valid in the vicinity of the 
instability point. If the eigenvalues of the corresponding unstable modes would differ significantly, the situation that all modes 
are in the unstable region would have the consequence that the maximum eigenvalue would be larger than zero, i.e. far away 
from the instability point. Thus, the order parameter equations would be no adequate approximation. Consequently, we only 
consider the case that Aaj,a]^ = J^x^jix^^ ^^T' >4^>4. ^ from which follows /J„ = , /f„ = VA^, A^, A^', A^' . 
However, it is possible that /Ju ^ . From now on the unstable modes are assumed to be given by 

A" = (1,0), (-1,0), (0,1), (0,-1). (72) 

This occurs, for instance, for monotonically decreasing cooperativity functions where we obtain the relation > /J (A 7^ 
0, ±1), by analogy with strings (see Sec. 1111 Ft . Then the maximum expansion coefficient is given by fx^ — /(a^^.a^') = fx^fxi^ 
for A" = 0, A2 — ±1 and A" = ±1, A2 = 0, respectively. Note, however, that the unstable modes ( |72| | could also arise for 
non-monotonic cooperativity functions as we will see below. 



C. Order Parameter Equations 

We specialize the order parameter equations (|9]l to planes. At first we determine the integrals (fTOl l. ( fTTT l of products of 
eigenfunctions ( |69] l, which read 

/ I \(«-l)/2 

-^A(i)A(2)... A(") ^ (Zli^j '5a(i)+A(2)+...+A("),A ) (73) 

f I \ ("-2)/2 

>^A(i)A(2)... A(") = ( j^^j^^ j '^A(i)+A(2) + ... + A("),0 • (74) 

For the order parameter equations ^ we need in (fT2] l. ( fTSl l. ( flST l the special cases 

-^A' A" = /T T '^A'+A",A , ^A' A" A'" = r T ^A' + A"+A"',A , J\' A" = ^A',-A" • (75) 
^ L1L2 1^11^2 

Note that the unstable modes jT^l i have the property A" + A" ^ A^. Thus, the quadratic coefficient ( fT2b vanishes due to the 
first integral of (|75]i: Alr\l,ll„ = 0. To yield a concise calculation of the order parameter equations, we introduce modes A" 
which are complementary to the unstable modes A" by permuting the two components: 

_ r (o,±i) if A" = (±1,0), 

~ \ (±1,0) if A" ^ (0,±1). 

The integrals (|75T l involve selection rules for the appearance of cubic terms (fT3T l which are analogous to those for strings. This 
condition turns out to be A" + A" + A" = A" for (|72] |. which leads to the following nine possibilities: 

(A"', A"", A"'") = (A", A", -A"), (A", -A", A"), (-A", A", A"), (A", A", -A"), (A", A", -A"), 

(A", -A", A"), (A", -A", A"), (-A", A", A"), (-A", A", A") . (77) 

In this way it can be shown that in total 14 possible cubic terms have to be taken into account. The resulting order parameter 
equations read 

Ux^x^^ = ci{Ux^xifU-x^-\i + C2Uxj^x^U-xi^\iUx-^-xi + csUxj^x^ (Ux^xiU_xj^_x^^ 
+U-x^-x-^U_x:^xi) +CiUx^xiUx^^-x^^U_xux^^+c^Ux^x^i^ {Ux^x^^U_x^-x^^ +'^-a^a]^^^a^-a]^) 

+C6C/-A^A^(7A.A;i'^A^-A]^ +C7C/a^AS (c^A^ A^^ ^^-A^-AJ^ + f^A^-A^C^-A^A^) 

+ C^A"A" t^A"-A" f/-A"A" + C^-A"A" t^A"-A" C^A"A" ) ■ (78) 



The coefficients ci-cs are fully determined by the expansion coefficients fx of the cooperativity functions and the control 
parameter a, as is documented in Ref. [SJ. 
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D. Real Variables 



To investigate how the complex order parameters contribute to the one-to-one-retinotopy between the planes, we have to 
transform them to real variables. To this end we introduce the transformation 



Thus, with the help of the function 



A^A^,A^ Aj^ ,Aj, A^ V ' / '^T^R Xj, X^ Xj. Xj^ '^T^R A^ X^ A^ A^ 

^T^H A" A" A" A" '^T^R A" A" A" A" 



the complex order parameter equations ( fTSl ) are transformed to real ones as follows: 

ux-^xi = ux^^xi [ci (u^. + «A^Ai^) + C2 ("a^-aj^ + «a^-A'a) 

+C7 (4^A^ +^^!^A^ +^^-A^ +«!^_aO] +^6/lA^-Aii,AJA5^,A^-A5iK^ 



+ C8 



^A^A^,A^A^,A^A^ ("^i ^} ^^A?^ — A^jA^^A^^jA^^ — A^^"^' ^) 



^ ^A^ — A^.A^ — '^^'■^T"^H ^} ^A^ — A^.A^ — '^fi''^T"^H 



(79) 



(80) 



(81) 



The equations for the amplitudes vxj^xi^ have an identical structure as they are obtained from (|8T|) by exchanging the variables 
u and V. It can be shown that the order parameter dynamics dSTI ) is governed by a potential |8]. However, as the corresponding 
expression for the potential is lengthy, we will not discuss it here explicitly. Instead, we investigate different analytical cases 
which depend on the number of non-vanishing modes. In particular, we are interested in the emergence of retinotopical ordered 
projections between the planes. 



E. Retinotopic Projections: Non- Vanishing Modes 



We start with the assumption that only the amplitudes Uj, vj of one mode are different from zero. We consider the case j — 2, 
but the other cases yield analogous results. The unstable part (|7]) reads 



U(t,r) = C/io,-ioexp 



■o ( ^1 ''I 



;7_io,ioexp 



'.If - If 



which is equivalent to 



U{t, r) ~ U2 cos 



' IT ~ I? 



V2 sm 



The real order parameter equations dSTT l reduce to 



ci 



^"'if "If 



Cl 



tt2 = + — U2(W2 + ^2) 1 ''^2 = AW2 + -^W2('«2 + ^2) 



(82) 



(83) 



(84) 



We obtain constant phase-shift angles, which was already discussed in the case of strings (see Sec lIIIEl l. With ^ — y^u^+v^ 
it follows ^ = + (ci/4)^^. For the stationary case this leads to ^ = or ^ = -^Z— 4A/ci. We are only interested in the case 
^ 7^ 0. This case corresponds to a retinotopy between ri and ti, respectively. Thus, we have a retinotopic order only in one 
dimension and not in the whole plane. 

Now we examine the question, if two modes are able to generate a retinotopic state in the plane. As a typical example we 
consider the case that U2, V2 as well as ug, f s remain. Thus, the unstable part (|7J has the complex decomposition 



U{t,r) = t/io,-ioexp i2n 
+J7oi,o-iexp i2n 



h r2 

^2 L2 



[7-10,10 exp 
+ C/0-1,01 cxp 



i27r 
-?;27r 



U_ _ ri_\ 
L\ Lfj 

^2 r2 

^2" ^2 



(85) 
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which due to ( [79] l corresponds to the real decomposition 



U{t,r) — U2COS 2n 
+us cos 2tt 
The real order parameter equations dSTT l read 



ti_ 

T 



ri 



Li Lf 



■ V2 sm 



t2 



r2 



L2 L2 



vs sm 



2tt 
2tt 



_ti n_ 



^2 



^•2 
^2 



(86) 



U2 = 


A^ 




h vl) - 


C7. 2 




U2 


^2 = 


A^ 










V2 


Ms = 


A^ 


Cl/ 2 

- jK- 




C7. 2 




Us 


Us = 


A^ 


Cl/ 2 

- jK- 




C7. 2 

f - 







(87) 



Again we obtain constant phase-shift angles. With the amplitudes ^ = y^u^+v^ , rj = \/v^+v^ the following coupled 
equations result 



(88) 



We investigate under which conditions the two retinotopic modes coexist. If ^, 7^ 0, we obtain from £, = — the relation 
(C^ — ?y^)(ci — C7) = 0. As we should minimize the restrictions for the coefficients ci, C7, we have in general ci 7^ C7, so we 
conclude £, — rj. Inserting this result in dSSl ) for the stationary case leads to 



4A 



Ci + C7 



(89) 



As the amplitudes ^, 77 have to be real and A > 0, the coefficients ci, C7 have to fulfill the condition ci + C7 < . Furthermore, 
we require that the coexistence of both modes is stable. To this end we consider the potential 



which reproduces according to 



? = ' ^ = 



drj 



(90) 



(91) 
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FIG. 8: The potential i90l for two coexistent retinotopic modes. According to d89l l there is a maximum at ^ = = and a minimum at ^ = r;. 
The right plot shows the equipotential lines, which pronounces the extrema of V. 




FIG. 9: The contribution of the center manifold S{U) leads to a more distinct concentration of the connection weight w{t,r) ~ 1 + U {t, r) + 
S{U {t, r)) around t = r, as compared with the approximation w(t, r) « 1 + [/ {t, r). 

the amplitude equations dSSl ). A stable state corresponds to a minimum of the potential ( |90l l. which leads to the condition 
ci < C7. With the explicit form of the coefficients ci and cy derived in Ref. |8] these considerations lead to the result, that either 
^20,20 Qj. ^11,11 could vanish. If 7^0,20 ^ jj^g stability is guaranteed by -57^^'^^/3 < 7 < -77^^'^^, whereas 711^11 = 
demands 57^°'^°/3 < 7 < —720.20 example we consider the first case and assume a special form of the cooperativity 

functions dTOb with the expansion coefficients 

/^io = /o';i = /±io = /o±i-0.1, /^i,±i = 0.05, /«,±i = -0.1, (92) 

thus it follows 7 = 0.01 and ^"^^'"^^ = —0.005. However, although the cooperativity function of the tectum is monotonically 
decreasing, the cooperativity function of the retina is not monotonic in this case, as is illustrated in Fig.|2l The corresponding 
potential ( |90l ) is shown in Fig. [8] 

F. Center Manifold 

In the next step we analyze the influence of the center manifold S{U) for the special case ( |92] l. The connection weight 
can be represented as 'w{t, r) = 1 + U{t, r) + S{U{t, r)) (see Sec. Using the relations ( fT4l i and ( fTSl l the stable part is 
approximated in the second order We rewrite the result to real variables and eigenfunctions, where we set ' = — 2. 

In this way we compare U {t, r) with U{t, r) + S{U{t, r)) in Fig. |9] It is evident that the retinotopic projection gets sharper 
due to the contribution of the center manifold, i.e. the projection is maximal around the point t = r. This corresponds to the 
situation found in the case of linear strings; the contribution of the higher modes have the tendency to support the emergence of 
retinotopic order. 
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w{i, f) 




retina tectum 



FIG. 10: The cells of retina and tectum, which are assumed to be continuously distributed on unit spheres, are represented by their unit vectors 
f and t, respectively. The two cell sheets are connected by positive connection weights 'w{t, r). 

V. SPHERES 

In this Section we apply our general model Ol] again to projections between two-dimensional manifolds. Now, however, we 
consider manifolds with constant positive curvature. Typically, the retina represents approximately a hemisphere, whereas the 
tectum has an oval form |2, 8]. Thus, it is biologically reasonable to model both cell sheets by spherical manifolds. Without 
loss of generality we assume that the two cell sheets for retina and tectum are represented by the surfaces of two unit spheres, 
respectively. Thus, in our model, the corresponding continuously distributed cells are represented by unit vectors f and t. Every 
ordered pair {i, f) is connected by a positive connection weight w{i, r) as is illustrated in Fig.fTO] The generalized Haussler-von 
der Malsburg equations ^ for these connection weights are specified as follows 

w{i,f)=f{i,f,w)-'^^ jd^r f{i',f,w)-^^^ Jdnr,f{i,f',w), (93) 
where the total growth rate is defined by 

fii,f,w) ^a + w{i,f) Jdflt' Jdnr.'CT{i-i')cR{f-f')w{i',f'). (94) 

The integrations in (|93] | and ( |94| | are performed over all points t, f on the spheres, where dilt, dflr represent the differential solid 
angles of the corresponding unit spheres. Note that the factors Svr in Eq. (|93] l are twice the measure of the unit sphere. After 
discussing the linear stability analysis around the homogeneous solution of the generalized Haussler-von der Malsburg equations 
(l93l l in Subsec. rV Al we perform the nonlinear synergetic analysis in Subsec lVBI which yields the underlying order parameter 
equations in the vicinity of the bifurcation. As in the case of Euclidean manifolds, we show that they have no quadratic terms, 
represent a potential dynamics, and allow for retinotopic modes. In Subsec. IV Cl we include the influence of higher modes upon 
the connection weights, which leads to recursion relations for the corresponding amplitudes. If we restrict ourselves to special 
cooperativity functions, the resulting recursion relations can be solved analytically by using again the method of generating 
functions, which was already applied in the case of strings. As a result of our analysis we obtain a perfect one-to-one retinotopy 
if the global growth rate a is decreased to zero. 

A. Linear Analysis 

According to the general reasoning in Ref. lit] we start with fixing the metric on the manifolds and determine the eigenfunc- 
tions of the corresponding Laplace-Beltrami operator Afterwards, we expand the cooperativity functions with respect to these 
eigenfunctions and perform a linear analysis of the stationary uniform state. 

1. Laplace-Beltrami Operator 

For the time being we neglect the distinction between retina and tectum, because the following considerations are valid for 
both manifolds. Using spherical coordinates, we write the unit vector on the sphere as i = (sin d cos sin d sin Lp, cos d). The 
Laplace-Beltrami operator for the sphere has the well-known form 
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whose eigenfunctions are known to be given by spherical harmonics Yim (x), i.e. 

A^^^Yi,r,{x) = -l{l + l)Yi„^{x), (96) 
which form a complete orthonormal system on the unit sphere. 

2. Cooperativity Functions 

The argument of the cooperativity functions c(i • x' ) is the scalar product x-x' which takes values between — 1 and + 1 . There- 
fore the cooperativity functions can be expanded in terms of Legendre functions Pi(x ■ x'), which form a complete orthogonal 
system on this interval IIOI 7.221.1]: 

c{x.i') = Y.^-fiPi{i.x'). (97) 
Using the Legendre addition theorem ifisll we arrive, for each manifold, at the expansion ((3) 

oo L oc I 



L=OM=-L 1=0 m=-l 

Note that the normalization of the cooperativity functions and the orthonormality relations lead to the constraints = = \. 



3. Eigenvalues 

According to Sec. |lll a linear stability analysis around the stationary uniform state leads to the eigenvalue problem of the 
linear operator (HJl. It has the eigenfunctions f) = Yj^j^j{i)Yi^^{f) and the spectrum of eigenvalues (|6]l. By changing the 

uniform growth rate a in a suitable way, the real parts of some eigenvalues (|6]l become positive and the system can be driven to 
the neighborhood of an instability. Which eigenvalues (|6]l become unstable in general depends on the respective values of the 
given expansion coefficients /J, /;^. If we assume monotonically decreasing expansion coefficients f-[, //^, i.e. 

l = ./o^>/f >/2^> •••>0, l = ./o''>/f >/2''>--->0, (99) 

the maximum eigenvalue in (|6]l is given by Amax = A^™ = —a + fi f^- Thus, the instability occurs when the global growth 
rate reaches its critical value ac = fTfi^- ^^^^ instability point all nine modes with = (1, 1) and A/" = 0,±1, 

m" = 0, ±1 become unstable, where we have introduced the index u for the unstable modes. 



B. Nonlinear Analysis 



In this subsection we specialize the generic order parameter equations of Ref. Uj] to unit spheres. We observe that the quadratic 
term vanishes and derive selection rules for the appearance of cubic terms. Furthermore, we essentially simplify the calculation 
of the order parameter equations by taking into account the symmetry properties of the cubic terms. We show that the order 
parameter equations represent a potential dynamics, and determine the underlying potential. 



1. General Structure of Order Parameter Equations 

The expansion of the unstable modes (|7]i reads U{i,f) = Uff ™ ^im" )^im" ) ^i^^, correspondingly, the contribution of 
the stable modes dHJ is given by S{t, f) = SfJ^Y^j^lt )Yj^{f ) . Note that the latter summation is performed over all parameters 
(i, I) except for (L", /") = (1, 1), i.e. from now on the parameters (i, I) stand for the stable modes alone. The resulting order 
parameter equations (|9]l read 

^M^m" _ ^jjM"m"^j^m",m"'m"" jjM"'m^' jjM""m"" 

. r>m" ,m"'m""m""' i-rM^'m^" TTM^"m"" TTM^"'m""' /-mrw 
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With the integrals (W^ 



^Lmm...n^) 



7m(i)m(2'...m("' 



the quadratic coefficients (fTSI) read A 



d^X Yi*^ix) Y;(l) „(1) (i) y;(2) „(2) (£)••• (x) , 

d^x Yini „(i) (i) Yi(2) „(2) (i) • • • y,(.) (x) 



(101) 
(102) 



m .m m 

M^'.M^'M 



„„ = Z]'^ j"*^ *^ 7™]^ ™ , whereas the cubic coefficients ( fTSI l are 



m ,m m m 



'M""M""' 



1,1 1 
•Jl 1 



c 7m m 



1 fT fR f jM",Kr' 
8^/1 /l ^^,111 

I J r/-7" /--R I fT -pRl tM".M^'M jm",m"'m 
+ '{[Jl Jl +/l/lJ-'l,lL -'1,1; ^47^ 

A I rR\ rm",m"'m , r x X /'1 1 tM" .M^' m] \ 

[^ + Jl ) ^l.ll +0/0 Om00m"m"' [^ + Jl) ^1,1 L J 



7-m .771 m 777 



'1,1 1 1 



X0M"A/"' Jl 1 



1 



LI 



The center manifold coefficients 



Mm,M"m,^M"'m'^' 
LI 



read 



LI 



2A- A 



rM.AI" AI"' 7-771,777 "777"' 1 / jM'^M'^' jm.m^'m"" r 



+^11 i^ii Oia 



(103) 



(104) 



2. Integrals 

The order parameter equations contain the following integrals: JJ^ ,/™'["™ ^I^^^ , ^™i7i ™ ™ . For the first integral 

one easily obtains J{" ™ = (—1)™ (5m'. -m"- Integrals over three and four spherical harmonics can be calculated using the 
Wigner-Eckart-theorem, 



47r(2/ + 1) 



whence for V = = 1 it follows 



.7«,m7« = (7(1 0,1,01/, 0)C(l,m',l,m"U, m). (106) 

^47r(2/ + 1) 

As the Clebsch-Gordan coefficients C(Zi, 0, 12, 0) vanish if the sum li + I2 + I3 is odd ifTlll . we obtain " =0. 
Thus, the quadratic contribution to the order parameter equations dlOOb vanishes, by analogy with Euclidean manifolds. Non- 
vanishing integrals ( 11061 ) can only occur for I = and / = 2. Furthermore, the integrals /™]^™ ™ follow from 1"^"^ ™ = 
(_iy77 +777 j-m ,-mm jjjfggj-^jg Qygj- fQ^J■ spherical harmonics can also be calculated using the Wigner-Eckart-theorem, and 
the result is 



/3 = |/"-/"'|7773 = -i3 V I 3 + j 

xC[l ,m ,/ ,m |i3,m3)/; ^. (107) 

Specialyzing (1107b to / = = /" = = 1 and taking into account (1105b leads to J™];™^™ ™ cx 5m'+m"+m"' ,m- Thus, we 
obtain the selection rule that the nonvanishing integrals fulfill the condition m' 
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3. Order Parameter Equations 

To simplify the calculation of the cubic coefficients (1103b in the order parameter equations jlOOI l. we perform some basic 
considerations which lead to helpful symmetry properties. To this end we start with replacing m" by — m" which yields 

ij^j^^]^ =^1111 ■ Correspondmg symmetry relations can also be derived for the other terms in 

( |103t . Therefore, we conclude that the order parameter equation for [/-a^^"-™" is obtained from that of f/*^"™" by negating 
all indices Af" and m" with unchanged factors. Thus, instead of explicitly calculating nine order parameter equations, it is 
sufficient to restrict oneself determining the order parameter equations for U^^, C/^", U^^, and U^^. The remaining five order 
parameter equations follow instantaneously from those by applying the symmetry relations, as is further worked out in Ref. f^. 

To investigate how the complex order parameter equations contribute to the one-to-one retinotopy, we transform them to real 
variables according to 



=i({/0i + [/0-i)/2, U7 = ({/i"-C/-i")/2, Us = i(i7i" + C/-i")/2. (108) 



The resulting real order parameter equations turns out to follow according to 

_ dV{{u,}) 



dui 



(109) 



from the potential 



V{{uj}) = -| II "I - y"o - P2ul{ul + ul) - I32ul{ul^ + ul) - Pzul{ui + ul + ul + ul) 
j=o 

— \/2/34Uo(uiU5U7 + U2U5U8 + U2UeUj + u^uqUj — uiugus — u^u^u-j — UiU^u^ — u^UQUg) 

-/35{ul - ul){uiU3 + U2U4) - /35(m? - ul){uiU3 - U2U4) - 2/35U7U8(uiU4 + U2U3) 

-2(35U5Ue{u2U3 ~ U1U4) + ^[(3e{ul + Ug) + (3e{u^ + ul)]{ul +ul + ul + uf) 

P"! / 2 , 2\t 2 I 2\ o / 2 I 2\/ 2 I 2\ 
-y (Ui + M2)(W3 + W4) - /?3K + "ejK + 

[{ul + ulf + (ul + ulf] +^{ul+ ulf + |(u? + ujf . (1 10) 

The dependence of the coefficients (3i and Pi on the expansion coefficients /; and the control parameter a is found in Ref. 11] . 
Naturally, a complete analytical determination of all stationary states of the real order parameter equations (|109l l is impossible. 
However, we are able to demonstrate that certain stationary states admit for retinotopic modes. 



4. Special Case 

To this end we consider the special case where ui,U2, M5, Me, uj, us vanish. Then the order parameter equations ( I109l l with 
(II 10b for the non-vanishing amplitudes uq, U3, U4 reduce to 

iio = Auo + 2Piul + 2(33{ul + ul)uo , 
U3 = Au3 + 2P3ulu3 + (3s{ul + ul)u3 , 

U4 — AU4 + 2/33U0U4 + /38(m3 + u|)m4 . (Ill) 

Due to the relation U3/U3 — U4/U4 one obtains constant phase-shift angles, i.e. it holds U3 cx U4. Therefore, the system of three 
coupled differential equations can be reduced to two variables. To this end we introduce the new variable 

i^^ul + ul, (112) 

which leads to 

MO = Amo + 2l3iul + 2/33C'mo , C = + 2P3ul^ + Psf ■ (113) 
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The stationary solution, which corresponds to a coexistence of the two modes, is given by 

9 A „ A 



(114) 



where we used the relation (3s = l^i + Ps f^- Demanding real amplitudes uq, £, leads to the coexistence condition /^s + /Jg < 0. 
Furthermore, we require stability for this state. Therefore we consider the corresponding potential V{uo, £,), which can be read 
off from (EOll and (dnii: 



A 



4 ^ ■ 



(115) 



Stable states correspond to a minimum of V, which leads to the conditions 2/33 ~ /^s > , /Sa — /Js > 0. If all three 
inequalities are valid, both the uq- and the ^-mode coexist. If we set U4 = 0, without loss of generahty, the solution reads in 
complex variables according to ( II 12l i 



00 



A 



Thus, the unstable part, specified in Sec. IV B Tl is given by 



A 



u{i,f) 



A 



Ps + Ps 

Using the Legendre addition theorem reduces (I117l i to 



U{i,f) = J- 



A 



P3 + Pi 



■Pi{i-f) 



(116) 



(117) 



(118) 



with Pi{i ■ f) = i ■ f. Thus, the unstable part is minimal, if i and f are antiparallel, i.e. the distance of the corresponding points 
on the unit sphere is maximum. Decreasing the angle between t and f leads to increasing values of U{t, f), and the maximum 
occurs for parallel unit vectors. This justifies calling the mode ( II I8I 1 retinotopic. 



C. One-to-One Retinotopy 

Now we investigate whether the generalized Haussler-von der Malsburg equations (|93] | describe the emergence of a perfect 
one-to-one retinotopy between two spheres. To this end we follow the unpublished suggestions of Ref. 1 12] and treat systemati- 
cally the contribution of higher modes. Because the Legendre functions form a complete orthogonal system for functions defined 
on the interval [—1, +1], their products can always be written as linear combinations of Legendre functions. This motivates that 
the influence of higher modes upon the connection weights, which obey the generalized Haussler-von der Malsburg equations 
(|93] l, can be included by the ansatz 

00 

w{a) ^J2'^2l + l)ZiPi{a), (119) 
1=0 

where the amplitudes Zi are time dependent. 



1. Recursion Relations 



Inserting (II 19l l into the generaUzed Haussler-von der Malsburg equations i93i and performing the integrals over the respective 
unit spheres leads to 



^(2/ + l)ZiP(((7) 



a 



l-^(2/ + l)ZiPi(a) 

1=0 

^(2; + l)ZiPi{a)J2{2l' + l)Zi,fj;fP[Pv{a) - Z,\ 



(120) 
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The products of Legendre functions occuring in (11201 ) can be reduced to linear combinations of single Legendre functions 
according to the standard decomposition ifTol 8.915] 

I 

Pi{a)Pi,{a) ^Y.^^^''-''P'+''-^k{cj), l<l' (121) 

fc=0 

with the coefficients 

{2l' + 2l-Ak + l)av^kakai-k (2fc-l)!! 

(2;' + 2/-2fc+l)a,+,„, ' - k\ ■ ^^^^^ 

Thus, contributions to the polynomial Pjicr) only occur iff the relation k = (/ + /' — /)/2 is fulfilled. Furthermore, using the 
orthonormality relation of the polynomials yields the following recursion relation for the amplitudes Zf. 



£(2r + i)z. 



{2l + l)Zi = a[5i^o-{2l + l)Zi]-{2l + l)Zi{Zl + iflf^Zl) 

I I" 

I'^O l"^0 k^Q 

oo V 

- (2^" + ^)Zi"fv'fu' E Av.y'.k5k,iv+i"~i)/2 

l"=V + l k=Q 



(123) 



Note that Eq. (I123I I cannot be solved analytically for arbitrary expansion coefficients ff, fj^ of the cooperativity functions. 
Therefore, we restrict ourselves from now on to a special case. 

2. Special Cooperativity Functions 

For simplicity we assume for the cooperativity functions (|97] l that /o = 1 , /i 7^ , /; = for I 7^ 0, ±1. With this choice the 
recursion relation ( |123t for I = reduces to 

Zo--(a + Zo' + 37^i)(^o-l), (124) 
where we have used again the abbreviation 7 — fi f^- For Z 7^ 0, by taking into account ( I122l i. we obtain 

Z, = -{a + Zl + i^ZDZi + Z,Zi + ^Zi-,+^l + l)Zi+, ^j25) 

The long-time behavior of the system corresponds to its stationary states. They are determined by Zq = 1 from ( 1124b . whereas 
( 11251 ) leads to a nonlinear recursion relation for the amplitudes Zi with I 7^ 0. However, by introducing the variable 

a + i-lZ^juf 

" = — o 7 ( \ — ' (126) 

6"iZi(u) 

this nonlinear recursion relation can be formally transformed into the linear one 

{l + l)Zi+i{u)^{2l + l)uZi{u) + lZi^i{u), 1>1. (127) 

Thus, determining the stationary solution of the nonlinear recursion relation ( 11251 ) amounts to solving the Unear recursion relation 
( 11271 ) for Zi{u) in such a way that the self-consistency condition ( 11261 ) is fulfilled. 

3. Generating Function 

To determine the amplitudes Zi (u) we calculate their generating function 

00 

E{x,u) ^Y^^)^^ ' (128) 
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where we have the normalization 

E{0,u) = Zo{u) = I. (129) 

Multiplying both sides of ( I127l i with x' and summing over I > 1 leads to an inhomogeneous nonlinear partial differential 
equation of first order for the generating function: 

dE(x u) 

(x^ - 2ux + 1) ^(u- x)E{x, u) + Zi{u)-u. (130) 

ax 



Using the normalization condition ( |129t its solution is given by 



^/x'^ - 2ux + l + x-u 
1 + [Zi (u) — u\ In 



4. Decomposition 

We now determine the unknown amplitudes Zi{u). From the mathematical literature it is well-known that the recursion 
relation (1127b holds both for the Legendre functions of first kind Pi{u) and second kind Qi{u), respectively ifToll . Thus, we 
expect that the generating function (|131| l can be represented as a linear combination of the generating functions of the Legendre 
functions of both first and second kind, which are given by [10, 8 .92 1 ] and lilft, 8 .79 1 .2] ; 

oo ^ 

Ep{x,u) = y^Pi{u)x' ^ , , ^ (132) 

Va; — 2ux + 1 



\J x"^ ~ lux + 1 + M — a; 

, 

Eq{x,u) = y^Qi{u)x' = , . (133) 

\Jx^ — 2ux + 1 



Indeed, taking into account the expUcit form of the Legendre function of second kind for / = 11 1311 

QoH = ^ln^, (134) 
2 ii — 1 

the generating function ( I131I I decomposes according to 

E{x,u) = {l + [Zi{u) - u]Qo{u)}Ep{x,u) - [^i(u) - u]Eq{x,u) . (135) 
Inserting (11321 ). ( I133l l and performing a comparison with ( 1128b then yields the result 

Zi{u) = {1 + [Zi{u) - u]Qoiu)}Pi{u) - [Ziiu) - u]Qi{u) . (136) 

Thus, the amplitudes Zi (u) turn out to be linear combinations of Pi (u) and Qi (u). To fix the yet undetermined amplitude Zi (u) 
in the expansion coefficients of (11361) . we have to take into account the boundary condition that the sum in the ansatz dl 191 ) has 
to converge. 



5. Boundary Condition 

Because the Legendre functions Pi{<t) do not vanish with increasing /, we must require that Zi{u) vanishes in the limit 
I oo. The series of Legendre functions of first kind Pi{u) with fixed u > 1 diverges for I oo according to fio', 8.917]: 
Po{u) < Pi{u) < P2{u) < . . . < Pn{u) < . . . , u > 1. The Legendre functions of second kind Qi{u), however, converge 
to zero (see Fig.fTTTi. Thus, performing the limit Z — ^ oo in Eq. ( 1136b and using the explicit form [13] Qi{u) — uQo{u) — 1 it 
follows that Zi{u) is fixed according to Zi{u) — Qi{u) / Qq{u). With this we obtain that the result (1136b finally reads 

^'W'^TTTT' ^^^^^ 
Qo[u) 



which is not valid only for / ^ but also for I — due to ( |129b . 
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FIG. 11: The Legendre functions of first and second kind Pi{u) and Qiiu) for u> 1. We have Pi{l) — 1, whereas Qi(u) diverges for u [ 1. 
Important for the boundary condition of Zi{u) is the different behavior for increasing values of I: Pi(u) diverges, whereas Qi{u) converges 
to zero. 



6. Connection Weight 
Inserting ( |137| i into ( 11191 ) yields the following solution for the connection weight: 

Qoiu) ^ 

Using the identity ^ 8.791.1] and (fT34] i. we obtain for the connection weight 



w{a) 



In 



u + 1 
u-1 



(138) 



(139) 



Note that integrating ( |139t over the unit sphere leads to 47r, i.e. the total connection weight coincides with the measure M. 

On the other hand we have to take into account that the self-consistency condition ( 11261 ) yields an explicit relation between 
the variable u and the control parameter a. Indeed, we infer from ( |126t and ( |137t the following transcendental relation between 
a and u 



- = In 



u+1 
u- 1 



2 In 



u + 1 
u-1 



(140) 



which is depicted in Fig.fT2k. 



7. Limiting Cases 

The limiting value of (11401 ) for u ^ cx) is determined with the help of the expansion ifiol 1.513] 

k=l 

and reads lim a = j. Thus, we conclude that the case u — > oo corresponds to the instability point etc = fi" fi , which was 

u — > oo 

obtained from the linear stability analysis in Sec. IV Al Correspondingly, using again ( 11401 ), we observe that the connection 
weight ( |139t coincides in the limit u ^ oo with a uniform distribution: lim w{(t) = 1 . Another biological important special 

case is M I 1, where we obtain from the transcendental relation ( 1140b lima = . Furthermore, considering the limit u j 1 in 

nil 

( 11391 ) for cr 7^ u, we obtain w{(j) 0. On the other hand, integrating ( |139t for u I I over a yields 2. Therefore, we conclude 
that the connection weight (|139|) becomes in this Umit Dirac's delta function: 

limu;(cr) = 4(5(cr- 1) . (142) 

Thus, decreasing the control parameter a means that the projection between two spheres becomes sharper and sharper (see 
Fig. [T2b). A perfect one-to-one retinotopy is achieved for a = when the uniform and undifferentiated formation of new 
synapses onto the tectum is completely terminated. 
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U)((t) 



03 < a2 < ai 



1 





1 



It a) — 1 







+1 



b) 



FIG. 12: a) Relation ( |140t between the control parameter a and the variable it. b) The connection weight for different values of the control 
parameter a. For decreasing values of a the connection weight around a = +1 is growing. In the limiting case a ^ the connection weight 
w{a) becomes Dirac's delta function l ll42t . 



In this paper we have explicitly applied our generic model for the emergence of retinotopic projections between manifolds 
of different geometry to one- and two-dimensional Euclidean and spherical manifolds. By treating retina and tectum as strings 
we generalized the original approach of Haussler and von der Malsburg where both were modelled as one-dimensional discrete 
cell arrays. This change from discrete to continuous variables is physiologically reasonable because of the high cell density in 
vertebrate animals. By using a continuous instead of a discrete model we emphasized that we are not interested in the dynamics 
of the single cell but in the evolving global spatio-temporal patterns of the system. Furthermore, continuous variables are 
helpful to describe retinotopic projections between manifolds of different magnitudes as we have seen by the example of two 
strings of different lengths. In case of discrete cell arrays with different cell numbers it is not clear what a perfect retinotopy 
means, whereas in the continuous case a perfect one-to-one projection can be described without offending the bijectivity of the 
projections. Finally, as the one-dimensional string model could only serve as a simplistic approximation to the real biological 
situation, we have also investigated under which conditions retinotopic projections between planar networks of neurons arise. 
Obviously, this increase of the spatial dimension rendered the synergetic analysis so complicated that we were only able to 
treat physiologically interesting special cases. While for strings retinotopy was only possible for monotonically decreasing 
cooperativity functions, we found that this is no longer true for planes. 

Applying our generalized Haussler- von der Malsburg equations |1] to strings and to spheres, led to remarkably analogous 
results. Both for one-dimensional strings and for spheres we have furnished proof that our generalized Haussler-von der Mals- 
burg equations describe, indeed, the emergence of a perfect one-to-one retinotopy. Furthermore, we have shown in both cases 
that the underlying order parameter equations follow from a potential dynamics and do not contain quadratic terms. However, 
in contrast to strings, spherical manifolds represent a more adequate description for retina and tectum. Therefore, the spherical 
case represents an essential progress in the understanding of the ontogenetic development of neural connections between retina 
and tectum. 
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